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ABSTRACT 

Scenarios  are  an  important  tool  in  the  strategic  planning  process,  and  are 
increasingly  used  in  both  the  Defence  and  business  world.  This  paper 
describes  some  potentially  useful  scenario  analysis  methods  for  systematically 
selecting  and  developing  future  scenarios.  The  processes  of  each  method 
are  illustrated  with  small  examples.  We  also  demonstrate  a  single,  flexible 
approach  to  combining  these  methods  using  a  typical  Defence  strategic 
planning  problem.  Some  general  guidelines  to  consider  when  choosing  and 
using  an  appropriate  scenario  analysis  method  are  also  discussed. 
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Some  Methods  for  Scenario  Analysis  in  Defence  Strategic 

Planning 

Executive  Summary 

Scenarios  are  an  important  tool  in  the  strategic  planning  process,  and  are  increasingly 
used  in  both  the  Defence  and  business  world.  Scenario  analysis  has  emerged  as  a  tool 
for  strategic  planning  when  the  future  is  perceived  as  surrounded  by  a  high  degree  of 
uncertainty  and  complexity.  Scenario  analysis  techniques  characteristically  synthesise 
quantitative  and  qualitative  information,  constructing  multiple  scenarios  or  alternative 
portraits  of  the  future. 

Scenario  analysis  consists  of  the  three  basic  stages:  (1)  Problem  analysis  to  come  up  with 
an  exact  definition  for  the  problem  of  the  investigation;  (2)  System  analysis  to  identify 
relevant  external  influences  on  the  problem  to  be  investigated,  and  (3)  Synthesis  process  to 
examine  the  existing  interdependencies  between  the  influencing  factors  and  to  establish 
alternative  scenarios. 

Problem  analysis  helps  all  experts  and  stakeholders  gain  a  common  understanding  of  the 
problem  at  hand.  Based  on  this  consensus  the  problem  can  be  further  bounded  and 
structured.  System  analysis  expresses  the  problem  as  a  system  of  inter-related  dynamic 
components  (subsystems),  with  the  system  itself  linked  to  its  external  environment.  From 
every  subsystem,  a  number  of  representative  influencing  factors  relevant  to  the  problem 
are  then  identified.  Synthesis  process  establishes  a  logical  and  systematic  way  for  scanning 
the  range  of  possible  scenarios  and  for  selecting  main  scenarios  or  choosing  a  set  of 
scenarios  that  includes  all  plausible  futures. 

A  variety  of  creative  methods  such  as  brainstorming,  brainwriting,  round  table  dis¬ 
cussion,  and  the  Delphi  technique  can  be  employed  in  the  first  two  analysis  stages. 
There  are  two  basic  methodologies  for  implementing  the  second  and  especially  the  third 
stage  of  the  scenario  analysis:  (1)  Non-Bayesian  method  (e.g.  Morphological  Analysis, 
Battelle  approach.  Field  Anomaly  Relaxation)  and  (2)  Bayesian  method  (e.g.  Cross-Impact 
Analysis).  Some  extensions  based  on  both  classes  are  also  developed  (e.g.  Battelle 
approach  with  Cluster  Analysis,  Cross-Impact  Analysis  with  Integer  Programming). 

This  report  describes  these  scenario  analysis  methods  and  provides  a  possible  way  for 
combining  them  into  a  single  more  flexible  approach.  The  processes  and  mathematical 
formulation  of  each  method  is  presented  in  this  paper,  as  well  as  with  a  discussion  of 
the  pros  and  cons  of  employing  these  methods.  The  combined  approach  is  illustrated 
with  a  typical  example  and  numerical  experiment.  This  approach  enables  the  scenario 
development  process  to  start  with  relatively  simple  information,  gained  from  experts  and 
problem-owners,  and  through  a  rigorous  and  transparent  process  identifies  a  manageable 
set  of  representative  or  balanced  scenarios. 

To  carry  out  the  combined  approach,  information  about  the  mathematical  formulation 
of  each  method  has  been  used  to  generate  algorithms  for  developing  computer  decision 
support  tools.  These  decision  support  tools  automate  complex  calculations  and  enable 
users  to  use  and  combine  the  techniques  without  in  depth  knowledge  of  the  mathematical 
aspects.  The  paper  concludes  by  emphasising  some  general  points  to  consider  when 
choosing  and  using  an  appropriate  scenario  analysis  method  in  Defence  strategic 
planning. 
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1  Introduction 

Scenario  analysis  has  emerged  as  a  tool  for  strategic  planning  [7]  when  the  future  is 
perceived  as  surrounded  by  a  high  degree  of  uncertainty  and  complexity.  Scenario  anal¬ 
ysis  techniques  characteristically  synthesise  quantitative  and  qualitative  information, 
constructing  multiple  scenarios  or  alternative  portraits  of  the  future. 

Scenarios  are  an  important  tool  in  the  strategic  planning  process,  and  are  increasingly 
used  in  both  the  Defence  and  business  world.  This  report  describes  some  scenario 
analysis  methods  and  provides  a  possible  way  for  combining  them  into  a  single  more 
flexible  approach.  The  combined  approach  is  illustrated  with  a  typical  example  and 
numerical  experiment.  This  approach  enables  the  scenario  development  process  to  start 
with  relatively  simple  information  (from  experts  and  problem-owners)  and  through 
a  rigorous  and  transparent  process  identifies  a  manageable  set  of  representative  or 
balanced  scenarios. 

After  providing  an  overview  of  scenario  analysis  methodology  in  section  2,  the  next 
three  sections  of  the  report  review  specific  methods.  The  processes  and  mathematical 
formulation  of  each  method  is  presented,  as  well  as  with  a  discussion  of  the  pros  and  cons 
of  each  method.  This  information  has  been  used  to  generate  algorithms  for  developing 
computer  decision  support  tools  (Section  6,  see  also  [5]).  It  should  be  noted  that  the 
list  of  methods  is  not  exhaustive,  and  other  techniques  are  available  in  the  list  of  related 
publications  [7, 16,  32],  What  is  presented  here,  however,  is  a  cross-section  of  techniques 
that  provide  rigorous  and  transparent  methods  of  generating  scenarios  for  complex 
strategic  planning  problems.  In  Section  6,  we  present  a  six-step  approach  to  combining 
methods,  which  is  demonstrated  by  an  illustrative  example.  The  report  finally  concludes 
by  emphasising  some  general  points  to  consider  when  choosing  an  appropriate  scenario 
analysis  technique  in  Defence  strategic  planning. 


2  Scenario  Analysis  Methodology:  An  Overview 

In  this  report,  we  consider  a  scenario  analysis  consisting  of  the  three  basic  stages  (cf.  [1]): 

1.  Problem  analysis:  to  come  up  with  an  exact  definition  for  the  problem  of  the 
investigation. 

2.  System  analysis:  to  identify  relevant  external  influences  on  the  problem  investigated. 

3.  Synthesis  process:  to  examine  the  existing  interdependencies  between  the  influenc¬ 
ing  factors  and  to  establish  alternative  scenarios. 

The  problem  analysis  stage  helps  all  experts  and  problem-owners  gain  a  common 
understanding  of  the  problem  at  hand.  Based  on  this  consensus  the  problem  can  be 
further  bounded  and  structured.  The  system  analysis  expresses  the  problem  as  a  system 
of  inter-related  dynamic  components  (subsystems),  with  the  system  itself  linked  to  its 
external  environment.  From  every  subsystem,  a  number  of  representative  influencing 
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factors  relevant  to  the  problem  are  then  identified.  The  synthesis  process  establishes  a 
logical  and  systematic  way  for  scanning  the  range  of  possible  scenarios  and  for  selecting 
main  scenarios  or  a  balanced  mix  of  scenarios. 

A  variety  of  creative  methods  such  as  brainstorming,  brainwriting,  round  table  discus¬ 
sion,  and  the  Delphi  technique  [14]  can  be  employed  in  the  first  two  analysis  stages.  There 
are  two  basic  methodologies  for  implementing  the  second  and  especially  the  third  stage 
of  the  scenario  analysis: 

•  Non-Bayesian  method  (e.g.  Morphological  Analysis  (MA)  [35],  Battelle  approach  [32], 
Field  Anomaly  Relaxation  (FAR)  [2,  3,  20,  21]) 

•  Bayesian  method1  (e.g.  Cross-Impact  Analysis  using  a  system  of  equations  [24,  25], 
or  Goal  Programming  (GP)  [4]). 

Some  extensions  based  on  both  classes  are  also  developed  (Battelle  approach  with  Cluster 
Analysis  [1],  GP  with  Integer  Programming  [12],  or  Fuzzy  Logic  [26,  27,  33]).  These  are 
described  in  detail  in  the  following  three  sections. 


3  Non-Bayesian  Method 

3.1  Morphological  Analysis  (MA) 

3.1.1  Description 

Morphological  analysis  was  first  developed  by  Fritz  Zwicky,  the  Swiss  astrophysicist  and 
aerospace  scientist  based  at  the  California  Institute  of  Technology,  during  the  Second 
World  War  [35].  Essentially,  morphological  analysis  is  a  method  for  identifying  and 
investigating  the  total  set  of  possible  relationships  contained  in  a  multidimensional 
problem  [23]. 

Firstly,  the  system  or  function  under  examination  is  broken  down  into  subsystems 
(components  or  dimensions).  In  this  breakdown  of  the  system,  the  choice  of  components 
is  critical  and  requires  considerable  thought  which  can  be  based  on  results  of  the  problem 
analysis  (e.g.  from  brainstorming).  The  aggregation  of  components  must  also  represent 
the  whole  system.  Too  many  components  avoids  a  clear  analysis;  conversely,  too  few 
makes  for  an  oversimplified  analysis.  Obviously  a  workable  compromise  must  be  found. 
It  should  be  noted  that  too  many  and  too  few  are  subjective  based  on  the  needs  of  the 
analysis. 

Each  component  can  take  on  several  configurations  (states  or  hypotheses).  A  given 
scenario  is  characterised  by  the  choice  of  a  specific  configuration  for  each  of  the 
components.  There  are  as  many  possible  scenarios  as  there  are  combinations  of 
configurations.  All  these  combinations  represent  the  field  of  possibilities,  called  the 
morphological  space. 

Whe  central  idea  of  the  Bayesian  method  is  to  elicit  the  likelihood  distribution  for  future  scenarios  to  be 
projected  from  the  experts  in  the  field. 
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Figure  1:  Scenario-building  through  Morphological  Analysis 


For  instance,  a  generic  example  of  scenario  building  through  morphological  analysis  is 
shown  in  Figure  1.  The  morphological  space  presented  is  composed  of  4  components; 
each  having  between  two  to  four  configurations.  This  enables  one  to  identify  a  number 
of  possible  combinations  (i.e.  96),  that  is  the  product  of  the  number  of  configurations 
(4  x  3  x  2  x  4). 

However,  certain  combinations  and  even  certain  families  of  combinations  are  unfeasible 
(e.g.  incompatibility  between  configurations).  The  next  phase,  therefore,  consists  of 
reducing  the  initial  morphological  space  to  a  useful  subspace,  by  introducing  exclusion 
factors  or  selection  of  criteria  (e.g.  economic,  technical),  from  which  the  relevant 
combinations  can  be  examined  (for  more  details  see  Section  3.2  or  Section  3.3). 


3.1.2  Application  Issues 

Although  the  method  has  been  used  in  technological  forecasting  [16,  19],  it  lends  itself 
well  to  the  construction  of  scenarios,  in  which  the  Social,  Technology,  Economics,  Ecology, 
Politics,  and  Values2  (STEEPV  model  -  the  six  themes  for  thinking  about  the  future  [15]) 
dimensions  can  be  characterised  by  a  certain  number  of  possible  states.  A  scenario  thus 
becomes  nothing  more  than  a  route,  a  combination  bringing  together  a  configuration 
for  each  component.  Morphological  analysis  therefore  stimulates  the  imagination  and 
enables  one  to  scan  the  field  of  possibilities  systematically. 

2The  selection  of  these  themes  will  depend  on  the  problem  being  solved. 
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The  first  limitation  of  MA  stems  from  the  choice  of  components.  By  leaving  out  a 
component  or  simply  a  configuration  that  is  essential  for  the  future,  one  runs  the  risk 
of  leaving  out  one  complete  facet  from  the  range  of  possibilities. 

The  second  limitation,  of  course,  stems  from  the  sheer  bulk  of  combinations  which  can 
rapidly  submerge  the  user.  Therefore,  MA  without  computer  support  severely  limits  the 
number  and  range  of  parameters  that  can  be  employed. 

In  the  next  two  sections,  we  will  investigate  two  forms  of  MA.  The  first  one  [21]  was 
developed  by  Rhyne  and  later  [2,  3]  used  by  Coyle,  Crawshay,  McGlone  and  Sutton 
where  the  phrase  "field  anomaly  relaxation"  is  applied  to  a  systematic  approach  for 
reducing  the  morphological  space.  Another  form  of  MA  mostly  appears  in  German 
literature,  developed  by  Battelle  Institute  in  Frankfurt  [32], 


3.2  Field  Anomaly  Relaxation  (FAR)  Analysis 

3.2.1  Description 

Field  Anomaly  Relaxation  (FAR)  [20]  uses  the  same  structure  as  MA  but  changes  its 
terminology.  A  problem  space  in  FAR  is  divided  into  various  fields  (called  'Sectors'), 
and  within  these  fields,  there  are  different  descriptions  of  possible  states  (called  'Factors') 
that  this  field  can  take.  A  configuration  is  formed  when  one  factor  from  each  sector  is 
combined,  forming  a  Whole  Field  descriptor  of  an  overall  condition  within  the  problem 
space.  A  part  of  the  process  for  FAR  is  giving  it  a  symbolic  name  (e.g.  the  OPTEC  model: 
the  letters  (symbols)  comprising  the  name  OPTEC  in  [30]  were  chosen  from  sector  names: 
political.  Physical,  Technical,  Economic  and  Cultural  environments)  and  then  using 
these  symbols  to  describe  the  configurations  (e.g.  O1P2T2E2C2  means  the  first  sector  takes 
the  first  factor  value  and  all  other  sectors  take  the  corresponding  second  factor  value). 

Some  configurations  (combinations  of  the  factors  of  the  different  sectors)  can  be  elimi¬ 
nated  as  being  not  plausible  to  occur  in  real  life  ( relaxation  of  anomalies)  which  reduces  the 
span  of  possibilities  of  future  developments.  This  process  is  called  filtering  in  FAR. 

The  first  filter  is  to  create  a  matrix  of  pairs.  Using  a  simple  scoring  yes/ no,  experts  answer 
the  same  question  for  each:  'Can  we  think  of  a  pattern  within  which  these  two  factors 
might  coexist?' 

When  a  matrix  is  fully  scored,  all  configurations  containing  inconsistent  pairs  (having 
pair(s)  scored  'no')  can  be  removed  from  consideration. 

The  second  filter  is  to  judge  the  wholeness  of  entire  configurations.  The  coherence  of  each 
whole  pattern  is  considered  as  certain  combinations  of  pairs  may  be  reasonable,  but  when 
put  together,  produce  a  scenario  which  is  not  sensible. 
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The  last  process  in  FAR  is  to  compose  scenarios.  The  surviving  whole  field  configurations 
are  positioned  on  a  tree,  in  which  nodes  represent  possible  future  states  and  branches 
represent  transitions  from  one  configuration  to  the  next.  Narratives  are  developed 
around  these  transitions  and  states,  which  then  are  formed  into  scenarios. 

A  scenario /future /'Faustian'  tree  is  developed  indicating  possible  pathways  from  con¬ 
figuration  to  configuration  and  endpoints,  ranging  from  positive  transitions  to  the  right 
side  of  the  graph,  and  negative  transitions,  on  the  left  hand  side  of  the  graph.  While  the 
vertical  dimension  represents  time  starting  from  present  and  spanning  up  to  considered 
future  times.  It  should  be  noted  that  neither  dimensions  represents  a  linear  scale. 

If  a  snapshot  of  a  possible  future  is  required,  then  a  single  configuration  may  be  expanded 
out  into  a  scenario.  If  a  potential  transition  from  the  current  state  to  a  future  state  is 
important  (possibly  to  investigate  shaping  strategies),  then  a  chain  of  configurations  may 
be  used  to  develop  the  scenario. 

Note  that  the  tree  presentation  is  intuitive  in  nature  and  a  matter  of  imaginary  (we  don't 
really  know  what  we  have  until  we  try!).  Early  designation  of  scenario  themes  may  ease 
scenario  composition  and  help  to  match  configurations  with  postulated  scenario  lines. 

Table  1:  Matrix  of  Pairs 


Ai  A2 

Bi  B2 

Ci  C2 

Di  D2  D3 

Ai 

a2 

Bi 

b2 

X 

G 

c2 

X 

D, 

d2 

d3 

X 

X 

Table  2:  Filtering  Inconsistent /Implausible  Configurations 


No. 

Configuration 

Inconsistent? 

No. 

Configuration 

Inconsistent? 

1 

AiBiCiDi 

13 

'  MClDi 

Yes 

2 

AiBiCiD2 

14 

\  ■  CiD2 

Yes 

3 

MBiCiM 

Yes 

15 

'■  HC,D3 

Yes 

4 

Yes 

16 

'  Bod, 

Yes 

5 

AiBHd2 

Yes 

17 

■  Bc2d2 

Yes 

6 

Yes 

18 

i  B  c2d3 

Yes 

7 

AiB2CiDi 

19 

A2B2CiDi 

8 

A2B2CiD2 

20 

A2B2CiD2 

9 

MbiCM 

Yes 

21 

A2B2CiD3 

10 

A2B2C2Di 

22 

A2B2C2Di 

11 

AiB2C2D2 

23 

a2b2c2d2 

12 

Yes 

24 

A^m ' 

Yes 

Let  us  illustrate  the  method  with  a  generic  simple  example.  Its  symbolic  name  is  ABCD 


5 


DSTO-TR-2242 


2035 


Figure  2:  A  Typical  Scenario  Tree 


(4  sectors).  Each  sector  consists  of  2  factors  except  the  last  one  (D)  having  3  factors.  There 
is  thus  23  x  3  =  24  possible  configurations.  Suppose  the  matrix  of  pairs  is  given  in  Table  1 
with  inconsistent  pairs  marked  with  a  'X'.  Table  2  lists  all  configurations.  Those  without 
inconsistent  pairs  (11  configurations  in  this  example)  can  be  investigated  as  a  whole,  to 
judge  the  coherence  of  the  complete  configuration. 

Assume  that  only  10  configurations  survive  after  the  second  filtering.  To  construct  a 
scenario  tree,  one  may  first  choose  a  configuration  for  the  current  state,  then  affix  each 
remaining  configuration  with  a  scenario  theme  and  place  them  on  the  graph  based  on  its 
timeframe.  A  typical  scenario  tree  may  be  arranged  as  shown  in  Figure  2  (for  illustrative 
purpose  only). 


3.2.2  Application  Issues 

As  FAR  is  a  form  of  MA,  all  pros  and  cons  of  MA  will  be  seen  in  FAR.  It  offers  a  systematic 
and  disciplined  approach  to  the  formulation  and  manipulation  of  data  and  ideas,  also  a 
way  of  reducing  the  number  of  scenarios  in  MA  to  a  manageable  set.  However,  there  is 
a  considerable  amount  of  time  and  effort  required  from  participants.  Several  workshops 
for  a  particular  FAR  analysis  is  a  common  practice  [28,  30].  Indeed,  FAR  inventor,  Russell 
Rhyne,  proposes  that  participants  should  spend  at  least  a  week  just  looking  and  thinking 
about  the  scenario  tree  [21,  22]. 
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The  filtering  process  is  a  tedious,  error-prone  task  when  done  by  hand.  A  computer 
support  tool  should  be  used  to  automate  some  of  FARs  processes.  Although,  analysing 
and  interpreting  data  and  results  must  be  cautiously  scrutinised  by  experts. 

Some  FAR  analysis  [2]  uses  a  very  conservative  approach  and  eliminates  combination 
pairs  only  if  there  is  consensus  within  the  groups  of  experts.  After  eliminating  all 
scenarios  that  contain  factor  value  pairs  that  are  declared  impossible,  the  process  is 
usually  terminated  (i.e.  only  first  pass  filtering). 

Lack  of  certainty  about  the  implausibility  of  particular  factor  value  pairs  will  cause  such 
pairs,  and  futures  containing  them,  to  be  retained.  Furthermore,  the  consistency  of 
each  whole  pattern  (carried  out  in  second  pass  filtering)  may  be  difficult  for  experts 
to  determine  with  confidence.  However,  it  may  be  possible  to  eliminate  them  when 
probability  estimates  are  applied.  This  will  be  explained  and  illustrated  in  the  Bayesian 
methods. 


3.3  Battelle  Approach 

3.3.1  Description 

The  'Battelle'  approach  from  Battelle  Institute  in  Frankfurt  [32]  defines  future  states  of  the 
environmental  subsystems  by  factor  outcomes.  Here,  the  outcomes  are  specified  mutually 
exclusively  but  exhaustively  for  certain  factors.  Its  structure  is  the  same  as  Morphological 
Analysis  (MA),  but  uses  different  terms.  A  comparison  of  main  terms  used  in  the  non- 
Bayesian  methods  is  listed  in  Table  3. 

Table  3:  Terms  Used  in  the  non-Bayesian  Methods 


Morphological  Analysis 

Field  Anomaly 
Relaxation  Analysis 

Battelle 

Approach 

Definition 

Subsystem/ Component/ 
Dimension 

Sector 

Factor 

an  environmental  subsystem 
representations  which  do  not 
have  one  clearly  defined  fu¬ 
ture  state 

State/ Configuration/ 
Hypothese 

Factor 

Outcome 

an  alternate  future  state  that 
subsystem  can  be  become 

Scenario 

Configuration 

Scenario 

a  combination  of  the  future 
state  of  the  subsystem 

As  a  non-Bayesian  method,  the  Battelle  approach  explicitly  does  not  use  probabilities; 
instead,  it  determines  the  interdependence  between  the  individual  outcomes  by  asking 
experts  to  rate  the  compatibility/plausibility  of  each  outcome  pair.  The  subjective 
estimates  are  called  compatibility  ratings,  which  are  expressed  on  a  scale  from  1  to  5.  A 
compatibility  rating  of  5  indicates  two  possible  occurrences  are  very  compatible,  and  a 
rating  of  1  indicates  they  are  not  likely  to  occur  together.  Values  of  2,  3,  and  4  represent 
increasing  compatibility. 
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The  output  is  a  range  of  compatible  scenarios  and  their  average  compatibility  values. 
One  may  choose  certain  scenarios  for  further  analysis  based  on  average  compatibility 
value  or  worst  compatibility  value  criteria. 

For  instance,  assume  the  scenario  in  a  future  environment  is  mainly  determined  by 
three  critical  factors  X\,x2  and  x3,  and  each  factor  consist  of  two  mutually  exclusive  but 
exhaustive  outcomes  Xy  and  xk/2,  k  =  1,2, 3. 

Table  4:  Factors,  Outcomes  and  Compatibility  Ratings 


Compatibility  rating 


Factor 

Outcome 

*1,1 

*1,2 

*2,1 

*2,2  *3,1  *3,2 

Xl 

*1,1 

- 

*1,2 

- 

- 

x2 

*2,1 

4 

3 

- 

*2,2 

2 

5 

- 

- 

x3 

*3,1 

3 

4 

2 

4 

*3,2 

3 

2 

5 

3 

Assume  further  that  the  subjective  estimates  (Table  4)  give  us  a  compatibility  rating 
(. ky ),  \/i,j  =  1, . . . ,  6  which  is  represented  as  a  triangular  matrix  ( ky  =  ky).  In  this  example, 
there  are  23  =  8  possible  scenarios  (outcome  combinations).  Their  worst  and  average 
value  compatibility  is  shown  in  Table  5. 

Table  5:  Scenarios  with  Worst/Average  Compatibility  Values 


Compatibility 


Scenario 

Set  Worst  value 

Average  value 

(1)  *1,1  *2,1  *3,1 

(4,3,2) 

2 

3.00 

(2)  xlrl  x2/2  x3/i 

(2,3,4) 

2 

3.00 

(3)  xu  X2,l  *3,2 

(4,3,5) 

3 

4.00 

(4)  xu  x2  2  x3/2 

(2,3,  3) 

2 

2.67 

(5)  x1/2  x2/1  x3/1 

(3,4,2) 

2 

3.00 

(6)  x1/2  x2/2  x3/2 

(5,2,  3) 

3 

3.33 

(7)  xi/2  x2/i  x3/2 

(3,2,5) 

2 

2.67 

(8)  x1/2  x2/2  x34 

(5,4,4) 

4 

4.33 

When  we  use  the  worst  compatibility  value  criteria,  the  top  three  scenarios  3,  6  and  8  are 
recommended  for  further  analysis  while  only  scenario  4  and  7  are  eliminated  from  the 
selection  process  if  using  the  average  compatibility  value  >  3  criteria. 


3.3.2  Application  Issues 

Like  FAR,  the  Battelle  approach  uses  MA  structure  to  break  down  problems  into 
factor  outcomes,  which  are  specified  mutually  exclusive  and  exhaustive.  Therefore, 
the  imcompatible/  inconsistent  outcome  combinations  (scenarios)  are  deleted,  and  those 
scenarios  with  high  compatibilities  remain. 
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The  Battelle  approach  uses  a  5  point  scale  to  obtain  the  compatibility  estimates  for  every 
possible  pair  of  factor  outcomes.  It  provides  a  more  flexible  scheme  for  experts  to 
rationalise  their  assessments,  but  may  also  take  longer  to  determine  their  ratings. 

As  a  non-Bayesian  method,  the  Battelle  approach  does  not  consider  the  probabilities  of 
outcomes,  therefore,  the  selected  scenarios  may  not  correspond  to  real  states  (scenarios 
may  have  very  small  probabilities  and  could  not  practically  be  a  basis  of  a  meaningful 
planning  effort).  We  will  consider  a  modified  approach  by  combining  all  techniques  later. 


4  Bayesian  Method  -  Cross-impact  analysis 

Much  of  the  writing  about  the  future  continues  to  be  in  the  area  of  speculation  [9], 
either  in  the  sense  of  constructing  scenarios  reflecting  the  writer's  expectations  of  likely 
developments  or  in  the  sense  of  advocating  changes  in  current  conditions  that  are 
thought  to  be  desirable.  But  in  the  long  run  the  message  of  such  speculative  endeavours 
remains  unconvincing  unless  there  is  a  parallel  development  of  analytical  methods  by 
which  to  derive  reliable  forecasts  and  consequent  estimates  of  the  relative  likelihood  of 
scenarios  of  the  future,  as  well  as  measures  for  attaining  desirable  futures  and  plans  for 
implementing  such  measures.  Knowledge  of  the  likelihoods  of  future  scenarios  is  needed 
for  planning  in  industry  and  government  [24],  Thus  the  central  idea  of  the  Bayesian 
method  is  to  elicit  the  likelihood  distribution  for  the  variable  to  be  projected  from  the 
experts  in  the  field. 


4.1  Model  Settings 

Cross-impact  analysis  techniques  [8,  9,  24,  25]  make  assumptions  about  the  future 
developments  in  the  environment  of  an  organisation,  in  which  certain  factors  influencing 
the  problem  under  investigation  are  identified  and  either  occur  or  do  not  occur3.  Starting 
with  a  finite  set  of  factors  x,  (i  =  1, . . . ,  n),  scenarios  consisting  of  combinations  of 
occurring  and  non-occurring  factors  are  constructed  in  the  synthesis  phases. 

Table  6:  List  of  Scenarios 


Scenario 

Probability 

Factor 

*2 

%n— 1 

Xn 

1 

yt 

1 

i 

i 

1 

2 

3/2 

1 

i 

i 

0 

N-  1 

Vn-\ 

0 

0 

0 

1 

N 

3/n 

0 

0 

0 

0 

3It  should  be  noticed  that  the  Bayesian  method  does  not  require  an  MA  structure  (i.e.  to  break  down 
problems  into  future  states,  which  are  specified  mutually  exclusive  and  exhaustive).  Modified  models  will 
be  presented  in  Section  5.1  to  make  the  Bayesian  method  applicable  to  any  problem  with  an  MA  structure. 
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The  combination  of  n  given  factors  X\, . ..  ,xn  produces  N  =  2"  scenarios.  Table  6  is  a 
scenario  list  where  the  likelihood  or  probability  of  a  scenario,  an  unknown,  is  denoted  by 
the  variable  ys,  s  =  1 ,N;  the  number  1  shows  that  the  factor  occurs  and  0  indicates 
that  it  does  not  occur. 

Two  most  recent  techniques  [4,  24,  25]  employing  cross-impact  analysis  have  modelled 
the  interdependencies  between  different  factors  in  the  form  of  conditional  probabilities 
p{i\j)  (the  probability  that  factor  xt  will  occur  given  that  factor  Xj  has  occurred). 
By  solving  systems  of  linear  equations  or  linear  programming  models,  the  scenario 
likelihoods  or  its  bounds  are  determined  from  the  marginal  and  conditional  probabilities. 
We  now  examine  the  mathematical  formulation  of  these  two  models  with  a  small 
numerical  example. 


4.2  Sarin's  Model  -  System  of  Equations 

In  Table  6,  let  an  be  the  column  vectors  of  0's  and  l's  ,  y  the  column  vector  of  the 

scenario  probabilities  ys  and  \f  the  corresponding  transposed  vector  of  y.  The  probability 
vector  y  must  satisfy  the  following  system  of  linear  equations  [24,  25]: 

yldi  =  p(i),  i  =  l,...,n 

yl (fli  a  (ij)  =  p(ij),  i  =  1, . . . ,  n;  j  >  i 

yl (al  A  a.j  A  ak)  =  p(ijk),  i  =  1, . . . ,  n;  j  >  z;  k  >  j  (1) 

yf(fli  A  (i2  A  . . .  f\an)  =  p(  12 . . .  n), 

N 

YjVs  =  1’  3/s  >  0,  Vs 

S=1 

where  the  'A'  operation  indicates  a  component  by  component  multiplication  of  two 
vectors. 

A  sequential  approach  to  determining  scenario  probabilities  ys  has  been  developed  as 
follows: 

•  Experts  firstly  provide  marginal  probabilities  p(i). 

Then  the  bounds  on  the  joint  probabilities  p(ij )  (where  p(ij )  =  p(i\j)p(j),  the 
probability  that  factors  x,  and  Xj  both  occur)  or  conditional  probabilities  p(i\j)  can 
be  calculated  using  standard  probability  theory. 

•  Experts  can  then  supply  additional  estimates  p(i\j)  within  these  particular  bounds. 

In  a  sequential  manner  the  bound  on  higher  ordered  joint  probabilities  are  com¬ 
puted  and  in  turn  used  by  the  experts. 
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Table  7:  Marginal  Probabilities 


Table  8:  Joint  Probabilities 


p(  0 

Value 

p(*;) 

Lower  bound 

Upper  bound 

p(  f) 

0.7 

P(  12) 

0.5 

0.7 

P(2) 

0.8 

p(13) 

0.0 

0.3 

P(3) 

0.3 

p(14) 

0.2 

0.5 

P(  4) 

0.5 

P(  23) 

0.1 

0.3 

P(  24) 

0.3 

0.5 

P(34) 

0.0 

0.3 

rank  the  scenario  in  order  of  their  likelihoods  and  select  probable /possible  scenarios  for 
further  analysis. 

For  example,  in  a  four-factor  scenario  analysis,  the  marginal  probabilities  of  occurrence 
for  each  factor  has  been  assessed  in  Table  7.  Using  the  following  conditions  on  the  joint 
probabilities  p(ij), 

max{0,  p{i)  +  p(j)  -  1}  <  p(ij )  <  min{p(z),  p(j)},  i  =  1,2, 3, 4  and  j  >  i, 
the  bounds  on  p(ij)  are  calculated  and  given  in  Table  8. 

Within  these  bounds.  Experts  can  now  specify  the  joint  probabilities  (e.g.  p(  12)  =  0.65, 
p(13)  =  0.25,  p(14)  =  0.25,  p( 23)  =  0.25,  p(24)  =  0.35  and  p(34)  =  0.1). 

Note  that  we  don't  illustrate  in  the  example  above  the  iterations  for  estimating  second- 
and  higher-order  joint  probabilities  as  the  next  model  using  Goal  Programming  (GP) 
formulation  will  remove  the  need  of  such  estimates  which  are  not  very  appealing  to 
participants. 


4.3  De  Kluyver's  &  Moskowitz's  Model  -  Goal  Programming 

This  model  only  requires  the  marginal  probabilities  p(z)  and  the  first-order  conditional 
probabilities  p(i\j)  with  its  bounds  p(z|/)~  and  p(i\j)+  (the  minimum  and  maximum 
of  p(z|;‘),  respectively).  Recognising  possible  inconsistency  between  the  estimates  of 
the  conditional  probabilities  and  the  marginal  probabilities,  the  model  using  Goal 
Programming  (GP)  formulation  seeks  to  find  consistent  values  from  which,  in  turn,  the 
probability  of  each  scenario  can  be  computed. 

It  introduces  the  theoretically  accurate  conditional  probabilities  p*(i\j)  which  fulfill  the 
axioms  of  probability  theory.  In  the  formulation,  the  inconsistency  is  written  in  the  form 

p*(*l;)  +&ij-6ij  =  viAi)' 

where  the  deviation  terms  Sjj  and  ib  measure  the  inconsistency  (the  difference  between 
the  theoretically  accurate  and  the  estimated  conditional  probabilities). 
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The  objective  of  this  GP  is  to  minimise  tl 
8  =f  max{(G,  }.  The  GP  is  of  the  form  [4]: 

minimise 
subject  to: 

l/tti 

y\cti  A  fly) 

N 

I> 

S=1 

p*(fl;)  +6ij  -< 

pO'l/T  <  p*01/) 

o  <  6;.  <  <5; 

ys  >  o, 

i  =  l,...,n;  j  > 


maximum  deviation  which  is  denoted  by 

8 

=  p(0/ 

=  p*(»I;)pO-)/ 

=  1,  (2) 

=  piAi), 

<  p(.Aj)+> 

0  <  8+j  <  8 

s  =  l 


Any  solution  to  GP  in  (2)  will  yield  a  set  of  the  scenario  probabilities  ys,  corrected 
conditional  probabilities  p*  (i\j)  and  their  deviations  8~,  8^  from  the  most  likely  assessed 
values  p{i\i)-  Upon  observing  the  results,  the  experts  may  wish  to  revise  the  assessment. 
This  method  is  illustrated  next  using  a  small  numerical  example. 


Suppose  that  the  four-factor  scenario  analysis  presented  in  Section  4.2  is  considered.  The 
marginal  probability  judgements  are  elicited  in  Table  7.  To  formulate  the  GP,  we  also 
need  to  estimate  the  most  likely  value  of  the  conditional  probabilities  p(i\j)  together  with 
its  bounds  and  p(i\j)+. 

The  probability  estimates  p(z),  p(j)  and  p(i\j)  to  be  consistent,  must  satisfy  the  following 
conditions 


max 


,  p( 0  +  p(j)  - 1 1 
p(j )  J 


<  p(AJ)  < 


min 


i  =  1, 2, 3, 4  and  j  >  i. 


From  the  above  condition,  one  can  thus  obtain  the  bounds  p(i\j)~  and  p(i\j)  1  ■  We 
suppose  further  that  p(i\j)  is  specified  within  these  bounds  and  given  in  the  last  column 
of  Table  9. 


Table  9:  Conditional  Probabilities 


Pairwise  factor 

(lj) 

Lower  bound 

P(Aj)~ 

Upper  bound 

P(Aj)  + 

Likely  value 

P(Aj) 

(T2) 

5/8 

7/8 

0.7 

(1,3) 

0 

1 

0.5 

(1,4) 

2/5 

1 

0.8 

(2,3) 

1/3 

1 

0.7 

(2,4) 

3/5 

1 

0.8 

(3,4) 

0 

3/5 

0.5 
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Using  any  mathematical  modelling  language  (e.g.  AMPL  [6]),  we  can  formulate  the  GP 
in  (2)  as  an  input  model  (see  Appendix  A  for  an  AMPL  implementation)  and  the  values  in 
Table  7  and  Table  9  as  input  data.  The  scenario  probabilities  ys  for  the  example  is  shown 
in  Table  10.  The  deviation  from  any  of  the  assessments  is  also  obtained  in  Table  11. 

The  Minmax  objective  (min  5)  is  used  in  the  first  case  which  results  the  change  in  all  most 
likely  values  of  the  conditional  probabilities  (third  column  in  Table  11)  and  the  likely 
scenarios  requiring  further  analysis  are  determined  as  1,  4,  7, 10  and  11. 


Table  10:  Scenarios  Probabilities 


Scenario 
(s)  X1X2X3X4 

Scenario  Probability  (ys) 
Case  1:  min  5 

Scenario  Probability  (ys) 

Case  2:  min  +  ^7  )  +  A4J 

(1)  1111 

0.135 

0.060 

(2)  1110 

0.000 

0.000 

(3)  1101 

0.040 

0.240 

(4)  1100 

0.325 

0.200 

(5)  1011 

0.053 

0.090 

(6)  1010 

0.000 

0.000 

(7)  1001 

0.110 

0.010 

(8)  1000 

0.038 

0.100 

(9)0111 

0.000 

0.100 

(10)0110 

0.113 

0.050 

(11)0101 

0.163 

0.000 

(12)  0100 

0.025 

0.150 

(13)  0011 

0.000 

0.000 

(14)  0010 

0.000 

0.000 

(15)  0001 

0.000 

0.000 

(16)  0000 

0.000 

0.000 

Table  11:  Conditional  Probabilities  &  Deviations 

Pairwise  factor 

{hj) 

Assessment 
P(*l  j) 

Case  1 

p*m 

Case  2 

p*m 

Case  1 

Case  2 

si 

Case  1 

Sb 

Case  2 

sii 

(L2) 

0.5 

0.625 

0.625 

0.125 

0.125 

0.000 

0 

(1,3) 

0.5 

0.625 

0.500 

0.125 

0.000 

0.000 

0 

(1,4) 

0.8 

0.675 

0.800 

0.000 

0.000 

0.125 

0 

(2,3) 

0.7 

0.825 

0.700 

0.125 

0.000 

0.000 

0 

(2,4) 

0.8 

0.675 

0.800 

0.000 

0.000 

0.125 

0 

(3,4) 

0.5 

0.375 

0.500 

0.000 

0.000 

0.125 

0 

The  second  solution  (Case  2)  reflects  the  expert 's  desire  to  keep  p*(i\j)  closer  to  the 
original  most  likely  estimates.  It  is  achieved  by  solving  the  GP  model  with  a  revised 
objective  function  +  S-- )  +  MS,  where  M  is  a  large  number,  say  10000.  So,  it 

minimises  first  the  maximum  deviation  5  from  any  of  the  most  likely  values  and  then 
secondly  the  sum  of  the  individual  deviations  Sh  and  .  We  can  see,  in  this  case,  the 
solution  (fourth  column  in  Table  11)  puts  only  p(l|2)  higher  than  its  estimate  (0.625 
versus  0.5)  and  the  likely  scenarios  are  changed  to  3,  4,  8,  9  and  12. 
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4.4  Application  Issues 

The  cross-impact  analysis  offers  a  way  for  modelling  and  assessing  the  interdependence 
among  several  relevant  factors.  The  difficult  simultaneous  consideration  of  all  relevant 
factors  can  be  avoided  and  sensitivity  analysis  for  each  factor  can  be  also  conducted. 

Cross-impact  analysis  requires  marginal  and  conditional  probabilities  for  the  pairs 
of  factors  as  input.  High  demands  are  therefore  placed  on  the  expert's  ability  and 
willingness  to  make  these  estimates.  The  process  of  assessing  conditional  probabilities 
for  factors  in  Sarin's  model  which  define  scenarios  is  both  complex  and  subtle.  Estimates 
from  experts  often  violate  probability  theory  axioms  and  consistency  tests  [18].  Although 
the  Goal  Programming  approach  can  guarantee  both  the  feasibility  and  consistency  of 
such  assessments. 

Cross-impact  analysis  outputs  a  ranking  of  scenarios  in  order  of  their  likelihoods.  This 
allows  us  to  select  possible  scenarios  for  further  analysis.  However,  these  methods  take 
all  scenarios  into  consideration.  In  consequence,  the  scenario  probabilities  are  often  very 
small.  It  is  suggested  that  a  cluster  analysis  (discussed  in  Section  5.2)  can  be  applied  to 
group  individual  scenarios  together  [17]. 


5  Extended  Methods 

5.1  Modified  Goal  Programming 

This  approach  [1]  allows  us  to  use  the  Goal  Programming  (GP)  formulation  (Section  4.3) 
for  problems  constructed  with  the  Battelle  approach  (Section  3.3). 

Experts  have  to  give  estimates  on: 

•  compatibility  ratings,  kjj  (as  in  Battelle  approach) 

•  marginal  probabilities  p(i)  on  the  occurrence  of  factor  outcomes  i  (as  in  Goal 
Programming  Approach;  however,  conditional  probabilities  for  the  pairs  are  not 
required  here!). 


5.1.1  Formulation 

Let  n  be  the  number  of  factor  outcomes  and  K  (where  K  <C  2")  be  the  number  of 
considered  scenarios,  the  modified  GP  is  of  the  form: 
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minimise 

YtfTj  +$)  +  MS 

kj 

(3a) 

subject  to: 

ylat 

<  p{  0/ 

(3b) 

yt{aj  A  fly) 

* 

VI 

(3c) 

K 

5=1 

<  i. 

(3d) 

=  p(y)/ 

(3e) 

p'W)  +  p'('J) 

=  p(  0/ 

(3f) 

VI 

1 

VI 

o 

0  <  5  J  <  5 

(3g) 

3/s  >  0, 

s  =  1  ,...,K 

(3h) 

i  =  1, ...  ,rr,  j  >  i  and  M  is  a  large  value,  say  10000, 

where  the  joint  probabilities  p(ij )  are  defined  by  the  transformation  of  the  marginal 
probabilities  p{i)  and  compatibility  values  kl]r  using  the  equations: 


p(v)  = 


p(i)p(j)  ~  V 

.p(0p(;')  +  ¥[''rP(')p(i)]' 


[hj-p(i)p(j)], 


if  1  <  k{j  <  3 
if  3  <  kij  <  5 


kj  <  p{ij)  <  Mi;/  i  =  1, .  •  ■ ,  n;  j  >  i 

def 

kj  =  max{0,  p{i)  +  p(j)  -  1}, 


uij  =f  min{p(i),  p(j)}. 


In  equation  (3e),  the  corrected  (or  final)  joint  probabilities  p*(ij )  of  the  preliminary 
(or  initial)  joint  probabilities  p(ij )  are  adjusted  by  deviation  variables  5  and  S  l  ;  5  is 
the  maximum  of  all  individual  deviation  variables;  and  p*(ij )  is  the  corrected  joint 
probability  that  outcome  i  will  occur  and  outcome  j  will  not. 

Some  other  notes  can  also  be  observed  on  the  modified  GP: 


•  The  objective  function  combines  two  functions  min^(<5;y  +  5^ )  and  minM<5  with 
the  latter  one  at  higher  priority  because  of  the  large  value  of  M. 

•  All  conditional  probabilities  are  replaced  by  their  correspondent  joint  probabilities 
so  that  constraint  (3c)  and  (3e)  are  basically  unchanged  if  comparing  to  GP  in  (2). 
Although,  only  a  subset  of  all  possible  scenarios  are  examined,  the  less-than-or- 
equal-to  constraints  are  used  on  (3b),  (3c)  and  (3d)  instead  of  equality  ones. 

•  In  equation  (3f),  the  probability  of  outcome  i  is  constrained  to  be  equal  to  the  sum 
of  the  joint  probabilities  for  outcome  i  and  every  other  outcome  both  occurring  and 
non-occurring. 
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5.1.2  Application  Issues 

The  modified  GP  model  provides  individual  scenario  probabilities,  but  because  of  the 
degenerate  solution  problem  in  linear  programming,  alternative  probabilities  exist.  It 
is  suggested  from  [1]  to  solve  the  modified  GP  first  to  obtain  the  minimum  possible 
deviation  (ni^ev)  and  then  to  create  a  new  objective  function  and  one  additional  constraint 
for  use  in  a  post-optimality  analysis.  Using  this  suggestion,  the  new  objective  function  is 


(4) 


Min  ys  or  Max  ys, 


and  the  additional  constraint  is 


YL^ij  +#ij)  +  MS  =  mdev,  Vz  =  1, . . . , tv,  j  >  i. 


(5) 


This  model  is  solved  for  each  of  the  K  scenarios  to  obtain  their  minimum  and  maximum 
probability  of  scenario.  The  arithmetic  mean  of  the  upper  and  lower  bound,  after  being 
adjusted  by  the  summation  of  all  scenarios  so  the  probabilities  summed  to  1,  defined  the 
probability  of  each  scenario. 

An  implementation  (in  AMPL)  of  the  modified  GP  formulation  (3a)-(3h)  together  with 
the  post-optimality  analysis  (4)— (5)  is  included  in  Appendix  B. 

5.2  Cluster  Analysis  -  Representative  Scenarios 

The  objective  of  scenario  analysis  is  to  develop  a  manageable  number  of  representative 
scenarios  that  can  be  used  in  strategic  planning.  The  optimal  number  of  scenario 
groupings  is  controlled  by  the  ability  of  the  end  user  (analysts,  experts,  stakeholders)  to 
conceptualise  the  alternatives  and  use  them  in  planning.  The  goal  of  finding  a  minimum 
number  of  scenarios  is  to  support  and  limit  the  work  of  the  scenario  writer  and  reader. 

The  idea  is  to  group  together  scenarios  that  are  'similar'.  This  raises  questions  of  how 
we  define  similarity  ( Measures  of  Similarity)  and  how  similar  do  they  have  to  be  to  be  put 
into  the  same  group  ( Cluster  methods).  Without  going  into  great  detail  of  cluster  analysis 
procedures  (see  e.g.  [13]  for  more  details),  we  will  only  examine  two  different  measures 
and  methods  which  have  been  suggested  and  used  in  the  strategic  planning  context  [1, 
17,  29]. 

5.2.1  Measures  of  Similarity 

The  definition  of  similar  is  subjective.  One  idea  is  to  use  the  distance  between  two 
scenarios. 

Squared  Euclidean  distance  for  use  with  the  scenarios  in  Bayesian  method  is  simply  the 
number  of  factors  on  which  two  scenarios  differ. 
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where  x?  and  xf,  i  =  1, ...  ,n  are  binary  value  of  the  factor  i  in  Scenario  p  and 
Scenario  cj  respectively  (see  Table  6  for  the  settings  of  Bayesian  method).  Note  that 
Manhattan  distance,  dist(s,  ^)  =  Y%=  1  |*f  —  *f  |,  gives  the  same  measure  in  this  case 
as  the  binary  setting  of  x1’  and  x‘j 

Compatibility  distance  for  use  with  the  scenarios  in  the  Battelle  approach  is  determined 
by  comparing  the  compatibility  ratings  between  the  factor  outcomes  in  one  scenario 
with  each  factor  outcome  in  another  scenario,  summing  all  of  these  compatibility 
levels,  and  dividing  by  the  number  of  factors  levels  compared.  The  resulting 
scenario  compatibilities  range  from  1  to  5. 

For  instance,  if  one  wants  to  calculate  the  compatibility  distance  between  Scenario 
(1)  *1,1X24X34  and  Scenario  (8)  *1,2X2,2X34  in  Table  5,  the  compatibility  values  of 
the  following  pairs:  (X14X24)  and  (*i,2*2,2);  (X14X34)  and  (*1,2X34);  (*2,1*34)  anc^ 
(*2,3*34)  must  be  compared  (absolute  value  of  their  difference),  then  summed  and 
divided  by  3.  By  taking  compatibility  values  from  Table  4,  we  obtain 

dist(Scenario  (1),  Scenario  (8))  =  — — ^  ^  3  ^  ^ ^  ~  1.33 


5.2.2  Cluster  Methods 

The  methods  used  in  scenario  analysis  fall  into  the  hierarchical  class  which  is  char¬ 
acterized  by  the  development  of  a  hierarchy  or  tree-like  structure  using  the  following 
algorithm: 

1.  Initialise  by  treating  each  scenario  as  a  separate  cluster  (with  only  one  member). 

2.  Compute  the  similarity  between  each  pair  of  clusters  ( cluster  similarity). 

3.  Find  the  two  most  similar  clusters  and  combine  them  into  a  new  cluster. 

4.  If  there  is  only  one  cluster  remaining,  stop.  Otherwise  go  to  step  2. 

To  run  this  algorithm,  we  need  to  define  similarities  between  clusters. 

Ward's  method  [34]  also  known  as  minimum  variance  method  (the  squared  Euclidean 
distance  to  the  center  mean),  aims  at  finding  compact,  spherical  clusters. 

Complete  Linkage  method  also  called  the  diameter  or  maximum  method  (based  on  the 
maximum  distance  between  scenarios,  one  from  each  cluster),  usually  performs 
quite  well  in  cases  when  the  scenarios  actually  form  naturally  distinct  'trends'. 


5.2.3  Application  Issues 

The  measures  of  similarity  and  the  cluster  methods  presented  here  are  not  the  only 
possibilities.  There  are  many  other  distances  (maximum,  Canberra,  manhattan,  binary, 
etc)  and  methods  (single,  average,  centroid,  etc)  in  Statistical  Modelling  Hl3l  that  can 
be  applied.  To  use  them  in  the  clustering  procedure,  we  should  ensure  that  they  are 
intuitively  satisfactory  to  the  user  (e.g.  the  ranking  of  the  assigned  distance  number  must 
agree  with  the  user's  judgment  of  the  relative  similarity  of  each  pair  of  scenarios). 
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Also  the  distance  measure  must  satisfy  a  set  of  metric  axioms: 

•  The  distance  from  a  scenario  to  itself  must  be  zero,  dist(p,  p )  =  0,  for  all  Scenario  p. 

•  The  distance  from  Scenario  p  to  Scenario  q  is  the  same  as  the  distance  from  Scenario 
q  to  Scenario  p,  dist(p,  q)  =  dist(^,  p)  for  all  p  and  q. 

•  For  any  three  scenarios  p,  q  and  r,  dist(p,  q)  +  dist(^,  r)  >  dist(p,  r). 

Most  statistical  packages  (e.g.  R  [10])  contain  a  cluster  analysis  module.  Once  a  distance 
and  method  are  selected  (or  defined  if  a  non-standard  one  is  used  such  as  compatibility 
distance  above),  the  clustering  procedure  can  be  carried  out  automatically  to  obtain  the 
cluster  tree  (called  a  cluster  dendrogram). 

Conceptual  or  practical  considerations  may  suggest  a  certain  number  of  clusters  (e.g. 
three  main  futures  possible  (Stable  and  Happy,  Economic  Recovery  or  Insurgency 
concern)  of  our  region)  or  the  theoretical  distances  (e.g.  at  least  a  chosen  distance  apart) 
at  which  clusters  are  combined  can  be  used  as  criteria. 

Interpreting  and  profiling  clusters  involves  examining  the  cluster  centroids  (mean  values 
of  the  cluster).  These  centroids  do  not  usually  correspond  entirely  to  possible  real 
scenarios  (e.g  'political  governance'  environment  (FAR  sector)  has  two  states  (FAR 
factor):  'stable'  with  value  1  and  'unstable'  with  a  value  of  2,  then  the  centroid  could 
be  given  the  value  1.7).  We  can  describe  the  scenario  with  a  vague  future  state  indicating 
the  range  of  variance  of  this  environment  in  this  cluster.  Alternatively,  mode  and  median 
of  the  state  within  a  cluster  are  examined  to  determine  a  representative  scenario  for  each 
cluster. 


5.3  Integer  Linear  Programming  -  Balanced  Mix  of  Scenarios 

Integer  Linear  Programming  (ILP)  formulations  [12]  are  developed  to  select  a  set  of 
scenarios  that  includes  all  future  states.  Selecting  a  minimum  number  of  plausible 
alternate  scenarios,  to  be  expanded  into  scenario  descriptions,  can  be  formulated  in  such 
a  way  that  each  state  (FAR  factor)  of  each  environment  (FAR  sector)  will  be  represented 
at  least  once  (or  twice,  or  three  times;  chosen  by  the  user). 

5.3.1  Formulation 

Denote  by  S,  the  set  of  all  scenarios  in  which  State  i  occurs.  Using  the  decision  variable  Z/c, 
taking  binary  value  0  or  1  according  to  whether  Scenario  k  (among  q  accepted  scenarios) 
is  selected  for  scenario  development,  the  ILP  can  be  written  as: 

<7 

Minimise  E  Zk 

k= 1 

subject  to  (6) 

Ez*>  Nir  Vi  =  1, . . . ,  n 

keSi 
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where  N,  is  an  integer  denoting  the  minimum  number  of  times  State  i  should  be  included 
in  a  scenario  definition. 

The  formulation  has  the  attraction  that  it  can  be  modified  and  extended  easily  by  adding 
a  variety  of  constraints  to  the  formulation.  For  example,  the  requirement  to  select: 

•  A  particular  scenario  can  be  represented  simply  by  setting  z/t  =  1  for  that  scenario. 

•  A  particular  combination  of  State  i\  and  State  h  to  be  at  least  R  times  could  be 
formulated  by  denoting  the  set  of  scenarios  that  contain  the  combination  S/lJ2  and 
adding  the  constraint  YLk eSt  t  zk  >  R- 

A  similar  formulation  results  if,  rather  than  requiring  a  state  to  be  represented  at  least 
Nj  times,  the  aim  is  that  the  total  probability  of  scenarios  in  which  State  i  occurs  is  set  to 
be  Pf.  This  obviously  requires  a  probability  estimate  for  scenario  k  as  input  data.  We 
can  use  the  arithmetic  mean  of  the  upper  (max  y^)  and  lower  (min  y*.)  values  probability 
estimates  in  the  Modified  Goal  Programming  (see  Application  Issues  in  Section  5.1)  for 
the  probability  Yj..  This  formulation  can  be  written  as: 

<7 

Minimise  E  zk 

k=  1 

subject  to  (^) 

J2ykzk>  Pi,  Vz  =  1, . . . ,  n, 

keSt 

where  Y^  =  ^(max  y ^  +  min  y^). 

Appendix  B  also  includes  an  implementation  of  both  ILP  (6)  and  ILP  (7). 


5.3.2  Application  Issues 

The  rationale  for  describing  each  particular  state  of  an  environment  in  a  scenario  context 
is  that  the  scenario  is  a  useful  way  of  exploring  the  significance  of  that  state.  At 
the  same  time,  it  is  pointless  to  describe  a  particular  state  in  the  context  of  other 
states  with  which  it  is  deemed  incompatible.  We  should  therefore  apply  the  technique 
after  selecting  compatible /plausible  scenarios  using  these  main  (non-Bayesian  or/and 
Bayesian)  methods. 

The  ILP  frequently  has  multiple  optima.  Alternative  solutions  should  be  found  and 
presented  to  the  end  user.  This  can  significantly  increase  the  flexibility  for  making  a 
decision.  For  finding  an  alternative  solution  of  an  ILP  problem  [31]  involving  only  binary 
variable  (z/c  6  {0, 1}  for  all  k),  we  just  add  the  following  constraint  to  exclude  an  existing 
solution: 

E2*  -  E  zk  ^  lBl  -  !/  (8) 

keB  keN 

where  B  =  {k\z^  =  1},  N  =  {k|z/c  =  0}  and  |B|  is  the  cardinality  of  set  B. 
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6  An  Approach  of  Combining  Methods  with 
Illustrative  Example  &  Decision  Support  Tool 

The  purpose  of  the  strategic  planning  process  is  to  reflect  possible  alternative  devel¬ 
opments  which  are  constructed  using  quantitative  data  as  well  as  the  experience  and 
intuition  of  Defence  experts  and  stakeholders.  However,  they  are  unlikely  to  be  interested 
in  the  mathematical  aspects  of  the  scenario  analysis.  Hence  information  required  from 
them  should  be  kept  as  simple  as  possible.  We  present  next  a  typical  example  in  Defence 
strategic  planning  and  an  approach  which  combines  all  the  above  methods  in  light  of 
these  requirements. 

As  seen  in  all  methods,  the  process  of  reducing  the  morphological  space  (non-Bayesian 
methods)  or  mathematical  models  (Baysian/Extended  methods)  may  appear  time  con¬ 
suming  and  complex,  a  decision  support  tool4  (DST)  is  also  designed  and  implemented 
using  the  familiar  spreadsheet  Microsoft  Excel  software  with  the  external  CPLEX  solver 
[11]  and  mathematical  modelling  language  AMPL  [6].  This  allows  all  numerical 
calculations  to  be  completely  automated. 


6.1  Typical  Example 

As  an  example  to  be  used  for  illustrating  the  combining  methods,  we  consider  the 
following  strategic  question  in  Defence  planning: 

Australia's  Joint  Operations  for  the  21st  century  states  regional  factors  (such  as  state 
fragility,  poor  governance  and  economic  underdevelopment)  may  affect  Australia's 
security  interests,  both  directly  and  indirectly.  As  a  result,  a  key  task  for  Australia's 
Defence  Force  is  to  contribute  to  a  stable  regional  environment. 

Contributing  to  a  stable  regional  environment  includes  being  able  to  defend  Australian 
territory  against  credible  threats  without  relying  on  the  combat  forces  of  other  countries, 
providing  joint  forces  to  contribute  to,  or  lead,  coalition  operations  in  Australia's 
neighbourhood  as  well  as  contributing  to  crisis  response  as  part  of  a  coalition  effort  in 
humanitarian  assistance  and  disaster  relief. 

This  leads  to  the  question,  what  will  Australia's  regional  environment  look  like  in  2030 
and  what  types  of  operations  will  Australia  be  required  to  respond  to  in  this  timeframe 
in  our  region? 


6.2  Six-step  Approach  of  Combining  Methods 

We  will  use  the  structure  of  the  non-Bayesian  methods  (Section  3)  to  break  down 
the  problem  space,  but  adopt  and  use  the  FAR  terminology,  'sectors'  and  'factors' 
(Section  3.2)  throughout  this  section.  Summary  of  the  approach  is  given  in  Table  12. 


4More  details  on  how  to  use  the  tool  including  a  full  package  (source  code,  test  files  and  documentation) 
can  be  obtained  by  contact  the  author. 
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Table  12:  Six-step  Approach  of  Combining  Methods 


Step 

Purpose 

Method 

1.  Description  of  future  states 

Identify  and  select  sectors  and  fac¬ 
tors  hypothesized  to  influence  the 
future. 

Brainstorm, 
MA,  FAR, 
Battelle 

2.  Assessment  of  states'  compatibilities 

Evaluate  compatibility /consistency 
values  between  pairwise  factors. 

Battelle 

3.  Determination  of  compatible  scenarios 

Define  a  criteria  for  plausible  /com¬ 
patible  scenarios  then  enumerate  all 
of  them. 

FAR, 

Battelle 

4.  Assessment  of  states'  possibilities 

Elicit  marginal  probabilities  on  the 
occurrence  of  factor. 

Bayesian 

5.  Analysis  of  scenarios'  possibilities 

Obtain  the  likelihoods  for  the  com¬ 
patible  scenarios  and  further  prune 
scenarios  due  to  their  likelihood. 

Modified  GP 

6.  Determination  of  main  scenarios 

Group  the  selected  scenarios  into  a 
few  main  ones  or  choose  a  balanced 
mix  of  plausible  futures. 

Cluster 
Analysis, 
1LP,  FAR 

Table  13:  Australia's  Regional  Environment  in  2030:  A  Morphological  Analysis 


Political 

Governance 

Pi:  Political  stability  in  most 
regions 

P2:  Unstable  political 
environment 

P3:  Collapse  or  change  in 
major  players 

Economic 

Growth 

Ei:  Developing 

E2:  Declining 

E3:  Collapse 

Social  Cohesion 

Si:  Tolerance  between 
groups 

S2:  Factionalisation  between 
groups 

S3:  Conflict  and  uprising 
between  group 

Implications  of 

S&T 

Ti:  Overwhelming  rate  of 
change  or  development  of 
technology 

T2:  Continuing  (comparable) 
advancement  of  technology 

T3:  Lagging  advancement 
of  technology. 

Health  and 
Habitat 

Hi:  Improving /Sustainable 

H2:  Degradation 

H3:  Collapse,  meltdown 

Type  of 
Operation 
required  by  ADF 

Ai:  Peacekeeping  /Peace 
enforcement  ADF  role 

A2:  Counter 
Insurgency/ Counter 
Terrorism 

A3:  Conventional  warfare 

A4: 

Humanitarian 

assistance 

ADF  Concurrent 
Obligations 

Ci:  Minor  commitment  to 
regional  Operations 

C2:  Major  commitment  to 
regional  Operations 

C3:  Commitment  to 
Operations  further  afield 

6.2.1  Description  of  Future  States 

The  first  step  in  developing  scenarios  is  to  identify  sectors  hypothesized  to  influence  the 
future  of  the  environmental  subsystems  investigated.  Although  the  number  of  sectors 
should  be  kept  to  a  minimum,  the  selected  sectors  need  to  be  comprehensive  enough  to 
reflect  all  relevant  concerns  about  the  future  and  be  thoroughly  defined  so  that  all  experts 
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understand  relevant  assumptions.  Six  to  seven  is  usually  recommended  for  the  number 
of  sectors. 

Two  to  five  possible  future  factors  are  designated  by  Subject  Matter  Experts  (SME)  for 
each  sector  using  historical  trends,  current  conditions,  and  expert  opinion.  These  factors 
are  mutually  exclusive  and  technically  exhaustive;  in  other  words,  other  factors  were 
thought  to  have  a  probability  of  occurrence  so  low  as  to  justify  their  exclusion. 

For  the  illustrative  example,  Australia's  Regional  Environment  in  2030,  the  description 
of  future  states  is  given  in  Table  13.  The  symbolic  name  is  chosen  as  PESTHAC  from  the 
7  sectors  (listed  in  the  far  left  boxes).  Each  sector  has  3  factors  except  Sector  A  (Type  of 
Operation  required  by  ADF)  which  has  4. 

6.2.2  Assessment  of  States'  Compatibilities 

The  interdependencies  between  factors  is  considered  in  this  step.  According  to  the 
Battelle  approach  (Section  3.3),  compatibility  ratings,  k,j,  are  given  by  asking  experts  to 
evaluate  how  compatible  two  factors  are. 

Table  14:  Sample  Data 


Compatibility 

P 

E 

S 

T 

H 

A 

C 

data 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

4 

1 

2 

3 

1 

0.3 

P 

2 

1 

0.2 

3 

1 

1 

0.5 

1 

5 

1 

3 

0.5 

E 

2 

1 

5 

3 

1 

0.4 

3 

1 

3 

2 

1 

1 

0.1 

1 

2 

3 

4 

3 

3 

2 

0.5 

S 

2 

4 

3 

2 

3 

3 

2 

1 

0.4 

3 

2 

1 

1 

1 

1 

1 

1 

1 

0.1 

1 

5 

3 

2 

3 

3 

1 

1 

4 

1 

0.6 

T 

2 

1 

3 

5 

3 

3 

1 

4 

2 

1 

1 

0.3 

3 

1 

2 

2 

2 

2 

1 

1 

1 

1 

1 

1 

0.1 

1 

2 

5 

3 

3 

3 

2 

5 

2 

2 

4 

3 

2 

0.8 

H 

2 

4 

2 

3 

3 

3 

2 

2 

5 

2 

2 

3 

2 

1 

0.1 

3 

1 

1 

2 

2 

2 

1 

1 

1 

1 

1 

2 

1 

1 

1 

0.1 

1 

2 

4 

3 

2 

4 

1 

5 

1 

1 

2 

4 

1 

5 

1 

1 

0.1 

2 

4 

2 

3 

3 

3 

1 

2 

4 

1 

3 

3 

1 

3 

3 

1 

1 

0.7 

3 

2 

2 

3 

1 

3 

1 

2 

1 

1 

2 

3 

1 

3 

1 

1 

1 

1 

0.1 

4 

1 

1 

2 

1 

2 

1 

1 

1 

1 

1 

2 

1 

2 

1 

1 

1 

1 

1 

0.1 

1 

2 

3 

2 

2 

4 

2 

2 

4 

2 

2 

3 

2 

2 

4 

1 

1 

5 

1 

1 

0.3 

C 

2 

3 

3 

4 

3 

3 

3 

3 

3 

3 

4 

3 

3 

3 

3 

3 

3 

3 

3 

3 

1 

0.5 

3 

4 

3 

3 

4 

2 

2 

4 

2 

2 

3 

3 

1 

4 

2 

1 

5 

1 

1 

1 

1 

1 

0.2 
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For  instance,  the  compatibility  data5  (below  the  diagonal)  used  to  illustrate  the  approach 
is  presented  in  Table  14  with: 

•  a  list  of  all  factors  i  (i  =  1,...,22)  corresponding  to  P\,  P2,  ■■■,C2  and  C3 

respectively 

•  the  compatibility  ratings  /c,;  for  every  two  factors  i  and  j,  where  kjj  E  {1, 2, 3, 4, 5}. 

Note  that  the  values  on  the  diagonal  (e.g.  p(E\)  =  0.5)  are  the  estimated  probabilities 
p(i )  for  the  individual  factor  i  which  will  be  described  in  Step  4  (Section  6.2.4). 

6.2.3  Determination  of  Compatible  Scenarios 

The  number  of  scenarios  are  exponentially  growing  with  the  number  of  factors.  Some 
combinations  of  factors  may  not  represent  plausible  scenarios.  In  order  to  decrease  the 
complexity  of  computation  and  consider  the  real  situations,  the  number  of  scenarios  are 
selected  by  the  following  rules: 

1.  A  compatibility  rating  between  any  two  factors  in  a  scenario  must  be  different  to  1 
(not  likely  to  occur  together),  and 

2.  The  average  of  individual  compatibilities  between  the  factors  in  each  scenario  is 
greater  than  or  equal  to  a  lower  limit  L,  or  the  number  of  compatibility  ratings  of  2 
(low  likelihood  of  occurring  together)  in  a  scenario  is  less  than  or  equal  to  an  upper 
limit  U,  where 

•  L  should  be  chosen  to  assure  the  remaining  scenarios  had  an  average  scenario 
compatibility  above  3  (in  other  words,  above  a  neutral  compatibility). 

•  U  should  be  below  half  the  number  of  the  sectors  in  a  scenario. 

Under  these  two  conditions,  scenarios  deemed  to  have  a  very  low  possibility  of  occurring 
are  eliminated.  In  some  cases,  the  participants  have  the  option  to  further  prune  to  a  subset 
of  these  compatible  scenarios  or  to  also  reintroduce  any  especially  interesting  scenarios 
which  were  excluded  due  to  their  incompatibility. 

In  the  illustrative  example,  all  compatible  scenarios  (i.e.  those  without  a  value  of  1)  are 
selected  using  the  following  value  of  L  and  U: 

•  Minimum  average  compatibility  value,  L  =  3.285. 

•  Maximum  number  of  “2  "  ratings,  U  =  3. 


5Note  that  all  data  presented  here  is  fictitious  and  used  for  illustrative  purposes  only. 
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Using  the  DST,  the  result  of  this  selection  process  is  shown  in  Table  15  where  the  'Factor' 
column  lists  all  accepted  scenarios  (e.g.  Scenario  6:  P2E2S1  T2H1  A|  C3). 


Table  15:  Selected  Compatible  Scenarios 


Scenario 

Factor 

Nunber  of 
Rating  2 

Avearge 

Compatibility 

Value 

Scenario 

Factor 

Nunber  of 

Rating  2 

Avearge 

Compatibility 

Value 

PESTHAC 

PESTHAC 

1 

1121121 

6 

3.286 

18 

3112113 

1 

3.714 

2 

1121122 

2 

3.381 

19 

3112122 

1 

3.286 

3 

1121221 

4 

3.524 

20 

3112222 

2 

3.143 

4 

1121222 

1 

3.524 

21 

3121122 

3 

3.095 

5 

2212112 

0 

3.667 

22 

3121222 

3 

3.143 

6 

2212113 

1 

3.810 

23 

3122122 

3 

3.048 

7 

2212121 

4 

3.286 

24 

3122222 

2 

3.190 

8 

2212122 

2 

3.238 

25 

3212112 

0 

3.619 

9 

2212132 

2 

3.238 

26 

3212113 

1 

3.714 

10 

2221  1  21 

4 

3.333 

27 

3212122 

1 

3.286 

11 

2221  1  22 

2 

3.286 

28 

3212132 

1 

3.286 

12 

2221221 

4 

3.333 

29 

3212222 

2 

3.143 

13 

2221222 

3 

3.190 

30 

3221122 

3 

3.095 

14 

2222122 

3 

3.095 

31 

3221222 

3 

3.143 

15 

2222221 

3 

3.333 

32 

3222122 

3 

3.048 

16 

2222222 

3 

3.095 

33 

3222221 

3 

3.333 

17 

3112112 

1 

3.524 

34 

3222222 

2 

3.190 

6.2.4  Assessment  of  States'  Possibilities 

This  approach  also  requires  marginal  probabilities  p(i)  on  the  occurrence  of  Factors  i. 
Because  possible  future  states  of  each  sector  are  considered  to  be  exhaustive  and  mutually 
exclusive,  the  assigned  marginal  probabilities  of  each  factors  in  each  sector  sum  to  1.  Also 
every  sector  usually  only  has  2  to  5  factors,  these  probabilities  are  quite  easy  to  elicit. 

The  marginal  probabilities  and  compatibility  ratings  obtained  above  are  then  used  to 
estimate  the  joint  probabilities  between  two  factors  and  to  serve  as  the  basis  to  obtain 
cross-impact  analysis  and  conduct  the  generation  of  scenarios. 

We  will  use  the  values  on  the  diagonal  of  Table  14  as  the  marginal  probabilities  for  the 
illustrative  example. 


6.2.5  Analysis  of  Scenarios'  Possibilities 
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We  now  calculate  the  probabilities  of  the  scenario  selected  in  the  previous  step  using 
the  modified  goal  programming  formulation  (5.1).  Here,  the  DST  will  call  the  external 
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CPLEX  solver  to  find  a  solution  for  the  modified  GP  (3a)-(3h)  which  is  modelled  through 
AMPL  language. 


Table  16:  Scenarios  Probabilities 


• 

eo 

-M 

'u 

•  i-H 

*Jh 

<75 

XI 

n 

c 

cS 

Q 

cs 

<L> 

rd 

0) 

x> 

u 

o 

u 

o 

CD 

$-1 

CD 

pH 

PH 

1 

5.90% 

10 

5.31% 

2 

8.55% 

11 

5.31% 

3 

3.24% 

12 

0.00% 

4 

2.95% 

13 

0.00% 

5 

2.95% 

14 

1.18% 

6 

2.36% 

15 

0.00% 

7 

3.54% 

16 

0.00% 

8 

3.54% 

17 

2.36% 

9 

1.77% 

18 

2.36% 

Scenario 

Probability 

Scenario 

Probability 

19 

8.85% 

28 

1.77% 

20 

0.00% 

29 

0.00% 

21 

6.49% 

30 

6.49% 

22 

2.95% 

31 

2.65% 

23 

1.18% 

32 

1.18% 

24 

1.18% 

33 

1.18% 

25 

2.95% 

34 

1.18% 

26 

2.36% 

27 

8.26% 

Based  on  the  solution  of  this  modified  GP,  the  upper  and  lower  bounds  for  all  selected 
scenario  probabilities  are  then  obtained  by  re-solving  the  modified  GP  with  new  objective 
functions  (4)  with  one  extra  constraint  (5).  The  arithmetic  mean  of  these  probabilities  is 
calculated  and  shown  in  Table  16. 

Note  that  Scenario  12,  13,  15,  16,  20  and  29  were  computed  to  have  probability  0 
throughout  the  parametric  analysis,  this  is  strong  indication  of  these  scenarios  being 
implausible.  So,  subject  to  expert  commentary,  these  scenarios  could  be  omitted  from 
further  consideration. 


6.2.6  Determination  of  Main  Scenarios 

In  the  final  step,  the  DST  is  used  again  to  perform  a  cluster  analysis  (Section  5.2)  for 
choosing  representative  scenarios  or  to  solve  ILP's  (Section  5.3)  for  a  balanced  mix  of 
plausible  scenarios. 

Table  17:  Sample  of  Cluster  Analysis 

Average 

Compability 

Cluster  Scenario  Value  Probability 


1,  4,  5,  7,  8,  9, 10, 

1 

14, 17, 18, 19,  21, 

22,  23,  24,  25,  27, 

28,  30,  31,  32,  34 

3.280 

77.0% 

2 

2,  26,  33 

3.480 

12.1% 

3 

3,  6, 11 

3.440 

10.9% 
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Cluster  Analysis 

Four  sets  of  trials  with  2,  3,  4  and  5  clusters  are  run  and  the  average  compatibility  rating 
for  all  scenarios  within  each  cluster  is  compared.  Table  17  shows  the  results  of  an  analysis 
with  three  clusters.  Table  18  displays  the  various  statistical  indicators  for  each  cluster  and 
also  proposes  a  representative  scenario. 

Note  that  the  basis  for  clustering  is  similarity  defined  by  a  distance  between  pairs 
of  scenarios  and  the  method  of  grouping  scenarios.  Here  we  use  the  inter-scenario 
compatibility  distance  and  the  complete  linkage  method  (Section  5.2)  as  to  maximise  the 
inter-scenario  compatibility  within  a  group  and  minimise  the  inter-scenario  compatibility 
between  groups.  However,  alternative  clustering  distance  and  method  are  possible. 

Table  18:  Cluster  Statistics  and  Representative  Scenarios 


Factor  Mean  Mode  Median  Maximum  Minimum  ^ePresenta^ve 

Scenario 


p 

2.43 

3 

3 

3 

1 

3 

E 

1.57 

2 

2 

2 

1 

2 

5— l 

S 

1.57 

2 

2 

2 

1 

2 

0> 

-M 

c n 

T 

1.64 

2 

2 

2 

1 

2 

u 

H 

1.29 

1 

1 

2 

1 

1 

A 

1.79 

2 

2 

3 

1 

2 

C 

1.86 

2 

2 

3 

1 

2 

P 

2.38 

2 

2 

3 

1 

2 

E 

1.77 

2 

2 

2 

1 

2 

(N 

J-H 

S 

1.54 

2 

2 

2 

1 

2 

-M 

c n 

T 

1.69 

2 

2 

2 

1 

2 

u 

H 

1.46 

1 

1 

2 

1 

1 

A 

2.00 

2 

2 

3 

1 

2 

C 

2.08 

2 

2 

3 

2 

2 

P 

2.75 

3 

3 

3 

2 

3 

E 

1.25 

1 

1 

2 

1 

1 

CO 

u 

S 

1.75 

2 

2 

2 

1 

2 

c n 

T 

1.75 

2 

2 

2 

1 

2 

3 

0 

H 

1.50 

2 

1.5 

2 

1 

2 

A 

1.75 

2 

2 

2 

1 

2 

C 

2.00 

2 

2 

3 

1 

2 

To  determine  whether  two,  three,  four  or  five  clusters  are  optimal,  an  average  compati¬ 
bility  rating  for  all  scenarios  within  each  cluster  is  calculated  to  determine  which  number 
obtains  a  maximum  average  compatibility  rating. 

The  representative  scenarios  (Table  18)  may  not  correspond  entirely  to  possible  real 
scenarios.  We  may  use  them  as  end-state  scenarios  and  others  in  their  cluster  as  transition 
scenarios  while  the  clusters  might  represent  different  branches  on  a  scenario  tree  (see  the 
last  step  of  FAR  analysis  in  Section  3.2). 
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Integer  Linear  Programming 


To  illustrate  another  possibility  for  selecting  a  minimum  number  of  plausible  alternate 
futures,  to  be  expanded  into  scenario  theme  descriptions,  we  use  the  selected  scenarios 
in  Table  15,  with  the  omission  of  implausible  Scenario  12,  13,  15,  16,  20  and  29. 

If  we  want  to  find  a  smallest  number  of  scenarios  that  cover  each  factor  twice  except  that 
the  factors  P3,  £3,  S3,  T3,  H3,  A3,  A4,  and  C3  (see  Table  13  for  the  description  of  the  factors) 
which  are  believed  to  be  insignificant  in  the  futures  (no  scenario  has  these  factors).  In  ILP 
model  (6),  Nt  thus  takes  the  value  {2,2, 0,2, 2,  0,2,2, 0,2,  2, 0,2,2, 0,  2,2, 0,0, 2, 2,0}  for 
each  factor  respectively.  The  DST  will  call  CPLEX  solver  to  find  a  solution  to  this  model. 
Scenario  3,  4,  6,  14, 18  and  33  are  listed  for  this  smallest  set  of  scenarios. 

If  the  total  probability  of  the  futures  in  which  the  factor  occurs  is  set  by  the  user,  then  ILP 
model  (7)  must  be  used.  For  example,  we  set  P,  (where  Y}--\  P  <  1)  to  the  value 


f—  —  0  —  —  0  —  —  0  —  —  0  —  —  0  —  —  00  —  —  01 

\  14'  14/  u/  14/  14/  14/  14/  14/  14/  14/  14/  u/  14/  14  /  u/  us  14  /  14  /  u  /  / 


respectively,  for  each  total  probability  of  factor  and  use  the  corresponding  scenario 
probabilities  listed  in  Table  16  as  the  values  for  Y/c.  The  DST  calls  CPLEX  solver  which 
outputs  Scenario  1,  3,  5,  6,  24  and  32  as  a  solution. 
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7  Conclusion 

There  is  probably  no  'best  method'  for  constructing  scenarios  in  Defence  strategic 
planning.  The  nature  of  the  problem  to  be  tackled  will  determine  the  most  appropriate 
method  to  be  used  in  that  instance.  However,  some  general  points  to  consider  are: 

•  Using  systematic  approaches,  strategic  planning  can  rationalise  the  often  ad  hoc 
process  of  selecting  futures  for  scenario  development. 

•  Morphological  analysis  stimulates  the  imagination  and  enables  one  to  scan  the  field 
of  possibilities  systematically. 

•  FAR  adds  to  MA  an  attractive  and  useful  way  of  formulating  and  manipulating  of 
data  and  ideas  by  presenting  a  scenario  or  'Faustian'  tree.  However,  this  requires  a 
considerable  amount  of  time  and  effort  from  participants. 

•  The  Battelle  method  is  straightforward,  simple  to  use  and  transparent. 

•  Plausibility  of  future  scenarios  can  be  captured  from  the  experts  in  the  fields  using 
the  Bayesian  method.  The  method  however  requires  intensive  information  input 
and  model  formulation  that  is  complex  and  must  be  implemented  in  a  special 
programming  language. 

•  An  attractive  feature  of  the  extended  and  combined  approach  is  the  fact  that 
starting  with  relatively  simple  information  from  experts  and  problem-owners,  main 
scenarios,  as  well  as  a  balanced  mix  of  plausible  futures  can  be  determined. 

•  A  computer  decision  support  tool,  similar  to  the  strategic  planning  tool  suite 
presented  in  this  report,  needs  to  be  used  to  automate  all  numerical  calculations 
in  each  step  of  the  combining  methods.  Although,  analysing  and  interpreting  data 
and  results  must  be  cautiously  scrutinised  by  experts. 
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Appendix  A  Modelling  Goal  Programming  with 

AMPL 

A.l  Model  File  "GP .  mod"  for  GP  Formulation 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 


set  N;  # number  of  factors 
set  NS;  ttnumber  of  scenarios 

set  NC;  ttnumber  of  first  —order  conditional  probabilities 

param  e{NS,  N}  binary ; 
param  aij{NS,NC}  binary ; 
param  p{N}  <=1,  >=0; 
param  pij{NC}  <=1,  >=0; 
param  paij{NC}  <=1,  >=0; 
param  pplus{NC}  <=1,  >=0; 
param  pminus{NC}  <=1,  >=0; 

var  y  {NS}  <=1,  >=0; 
var  ypij  {NC}  <=1,  >=0; 
var  dplus{NC}  <=1,  >=0; 
var  dminus{NC}  <=1,  >=0; 
var  d  <=1,  >=0; 

minimize  deviation_casel  :  d; 

minimize  deviation_case2  :  sum{  i  in  NCI  ( dplus  [  i  ] +  dminus  [  i  ])  +  10000*d; 

subject  to  Marginal.P  {j  in_  N}: 

sum{  i  in  NS}  y  [  i  ]  *  e  [  i  ,  j  ]  =p  [  j  ] ; 

subject  to  Sum  : 

sumj  i  hi  NS}  y  [  i  ]  =  1 ; 

subject  to  Conditional_P  {k  in_  NC}: 

sum{  i  in  NS}  y  [  i  ]*  aij  [  i  ,k]— paij  [k]*  ypij  [k]  =  0; 

subject  to  Conditional_Bound  {k  in  NC}: 
pminus[k]  <=  ypij  [k]  <=  pplus[k]; 

subject  to  Most_likely _P  {k  in_  NC}: 

ypij  [k]  — dplus  [k]  +  dminus  [k]=  pij  [k]; 

subject  to  def.dplus  {k  hi  NC}: 
d  —  dplus  [k]  >=  0 ; 

sub  j  ect  to  def_dminus  {k  hi  NC}: 
d  —  dminus  [  k  ]  >=  0 ; 

s _ > 
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A.2  Command  Script  File  "GP .  run"  for  Running  Model 


1 

2 
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4 

5 

6 

7 

8 
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10 

11 

12 

13 

14 
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16 
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19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 
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40 


model  GP .  mod ; 
option  show.stats  1; 

table  Tpii  "ODBC"  "testde.xls"  "pii": 

N  <—  [pii],  p  ~  value  IN; 

read  table  Tpii; 

table  Tp  "ODBC"  "testde.xls"  "p": 

NC  <—  [p],  pij  ~  pij  IN,  paij  ~paij  IN, 
pplus  ~pplus  IN,  pminus~pminus  IN; 

read  table  Tp ; 

table  Teij  IN  "ODBC"  "testde.xls"  "eij": 

{i  in  NS}  <-  [  eij  ] ,{  j  in  N}  <  e  [  i  ,  j  ]  ~  ( j  )  >; 

read  table  Teij  ; 

table  Taij  IN  "ODBC"  "testde.xls"  "  a  i  j  " : 

(i  in  NS}  <-  [  aij  ]  ,{  j  in  NC}  <  a i j  [i  ,  j  ]  ~  ( j  )  >; 

read  table  Taij; 

display  p,  pij  ,  paij  ,  pplus  , pminus  ,e  ,  aij  ; 

# - # 

#  Case  1:  Minimise  the  overall  maximum  deviation 

#  - # 

problem  Deviation_Casel  :  y,  ypij  ,  dplus ,  dminus ,  d, 
deviation_casel  ,  Marginal.P  ,  Sum,  Conditional _P 
ConditionaLBound  ,  Most_likely _P  ,  def_dplus  ,  def_dminus  ; 
solve ; 

display  y,ypij  ,  dplus  , dminus ; 

# - # 

#  Case  2:  Minimise  the  overall  maximum  deviation 

#  and  the  sum  of  individual  deviations 

#  - # 

problem  Deviation_Case2 :  y,  ypij,  dplus,  dminus,  d, 
deviation_case2  ,  Marginal.P,  Sum,  Conditional _P 
ConditionaLBound  ,  Most.likely  _P  ,  def.dplus  ,  def.dminus  ; 
solve ; 

display  y,ypij  ,  dplus  , dminus ; 


Note  that  the  values  in  Table  7  and  Table  9  are  stored  in  the  Microsoft  Excel  file 
"testde.xls"  to  generate  input  data  set  according  to  the  GP  formulation  described  in 
Section  4.3.  More  details  on  how  to  implement  Excel  macro  to  produce  the  input  data  set, 
including  a  full  package  (source  code,  test  files  and  documentation)  can  be  obtained  by 
contact  the  author. 
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Appendix  B  Modelling  Modified  GP  and  ILP  with 

AMPL 

B.l  Model  File  "mGP-MIP .  mod" 


set  N;  #  number  of  factors 
set  NS;  #number  of  scenarios 

set  NC;  #number  of  first-order  joint  probabilities 

set  NO;  ttnumber  of  outcomes  both  occuring  and  non—occuring 

param  e{NS,  N}  binary ; 
param  aij{NS,NC}  binary ; 
param  p{N}  <=1,  >=0; 
param  pijjNC}  <=1,  >=0; 
param  ppij{NO}  <=1,  >=0; 
param  apij{NO,NC}  <=1,  >=0; 
param  OF; 
param  ymax{NS}; 
param  ymin  {NS  } ; 
param  t ; 


var  y{NS}  <=1,  >=0; 
var  ypijjNC}  <=1,  >=0; 
var  dplusjNC}  <=1,  >=0; 
var  dminusjNC}  <=1,  >=0; 
var  d  <=1,  >=0; 


4 

4 


subject  to  MarginaLP  {j  hi  N}: 

sum{  i  in  NS}  y  [  i  ]  *  e  [  i  ,  j  ]  <=p  [  j  ] ; 

subject  to  Joint.P  {k  hi  NC}: 

sum{  i  hi  NS}  y  [  i  ]  *  a i j  [  i  ,  k]—  y p i j  [k]  <=0; 

subject  to  Sum  : 

sum!  i  in  NS}  y[i]<  =  l; 

subject  to  Most.likely  _P  {k  in  NC} : 

ypij  [k]  —  dplus  [k]  +  dminus  [k]=  pij  [ k  ] ; 

subject  to  P_occurring  {k  in.  NO}: 

sum{j  in  NC}  ap  ij  [k ,  j  ]  *  ypij  [  j  ]=  ppij  [k  ] ; 

subject  to  def.dplus  {k  hi  NC}: 
d  —  dplus [k]  >=  0; 


# - 

#  Modified  GP  formulation 

#  - 

minimize  deviation  : 

sum!  i  in.  NC}  (  dplus  [  i  ]  +  dminus  [  i  ])  +  10000*d; 
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sub j ect  to  def_dminus  {k  hi  NC}: 
d  —  dminus  [  k  ]  >=  0 ; 


# - # 

#  Post— optimality  analysis 

#  - # 

maximize  P_max{  i  hi  NS }  :  y  [  i  ] ; 
minimize  P_min {  i  hi  NS }  :  y  [  i  ] ; 

subject  to  deviation.extra  : 

sum{  i  in  NC}  ( dplus  [  i  ] +  dminus  [  i  ] )  +  10000*d=  OF; 


# - # 

#  Integer  Linear  Programming 

#  - # 

set  S { j  hi  N}:  = 

{k  iri  NS:  e[k,j]  =  l};  ttnumber  of  scenarios  in  which  outcome  j  occurs 

param  Nj{N}; 
param  Pj{N}; 
param  yave{NS}; 

var  z{NS}  binary ; 

minimize  No_Scenarios  :  sum  j  k  in  NS}  z[k]; 


subject  to  Factor_Value  (j  hi  N}: 
sum{  i  in  S[  j  ]}  z[i]>=Nj[  j  ]; 

subject  to  Fact  or  _P  {j  hi  N}: 

sum{k  hi  S[  j  ]}  yave  [k] * z  [k]  >=  Pj  [  j  ] ; 


B.2  Command  Script  File  "mGP-MIP .  run" 
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model  mGP-MIP. mod; 

option  show.stats  1; 
option  presolve  0; 

table  Tpi  "ODBC"  "test.xls"  "pi": 

N  <—  [pi],  p  ~  value  IN; 

read  table  Tpi; 

table  Tpij  "ODBC"  "test.xls"  "pij 
NC  <—  [pij],  pij  ~  value  IN; 

read  table  Tpij; 

table  Tppij  "ODBC"  "test.xls"  "ppij 
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NO  <—  [ppij],  ppij  ~  value  IN; 

read  table  Tppij ; 

table  Teij  IN  "ODBC"  "test.xls"  "eij": 

{i  in  NS}  <-  [  eij  ] ,{  j  in  N}  <  e  [  i  ,  j  ]  ~  ( j  )  >; 

read  table  Teij  ; 

table  Taij  IN  "ODBC"  "test.xls"  "aij": 

(i  in  NS}  <-  [  aij  ]  ,{  j  in  NC}  <  aij  [i  ,  j  ]  ~  ( j  )  >; 

read  table  Taij; 

table  Tapij  IN  "ODBC"  "test.xls"  "apij": 

(i  in  NO}  <-  [  apij  ]  ,{  j  in  NC}  <  apij  [i  ,  j  ]  ~  ( j  )  >; 

read  table  Tapij; 


# - # 

#  Solving  Modified  GP 

#  - # 

problem  Deviation_Min  :  y,  y p i j  ,  dplus  ,  dminus ,  d, 
deviation,  MarginalJP  ,  Joint_P  ,  Sum, 

Most_likely_P  ,  P_occurring  ,  def_dplus  ,  def_dminus  ; 

solve  Deviation_Min ; 


# - # 

#  Carrying  out  Post  — optimality  analysis 

#  - # 

problem  P_Max{i  in.  NS}:  y,  ypij  ,  dplus,  dminus,  d, 

P_max[i],  deviation.extra  ,  Marginal.P  ,  Joint_P  ,  Sum, 
Most_likely_P  ,  P_occurring  ,  def_dplus  ,  def_dminus  ; 

problem  P_Min{i  in.  NS}:  y,  ypij  ,  dplus,  dminus,  d, 

P_min  [  i  ] ,  deviation_extra  ,  Marginal.P  ,  Joint_P  ,  Sum, 
Most_likely_P  ,  P_occurring  ,  def_dplus  ,  def_dminus  ; 

let  OF:  =  deviation  ; 

display  OF,y; 

let  t :  =  0; 

for  {  i  in.  NS} { 

solve  P_Max [  i  ] ; 

let  ymax [  i  ]  :  =  y  [  i  ] ; 

solve  P_Min  [  i  ] ; 

let  ymin  [  i  ]  :  =  y  [  i  ] ; 

let  yave  [  i  ] :  =  ymax[  i  ]+ymin  [  i  ] ; 

let  t  :=  t+yave [ i ] ; 

}; 


display  {  i  hi  NS}  (yave  [  i  ]  /  t  ,  yave  [  i  ]  ,ymin  [  i  ]  ,ymax[  i  ] )  , 
sum!  i  in_  NS }  ( ymax  [  i  ]  +  ymin  [  i  ] )  /  2  , 
sum! i  hi  NS}  yave[i]/t; 
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# - # 

#  Integer  Linear  Programming  (outcome  represented  Nj  times) 

#  - # 

table  TF  "ODBC"  "test.xls"  "F": 

N  <—  [F],  Nj  ~  Nt  IN,  Pj  ~  Pr  IN; 
read  table  TF ; 
display  Nj ,  Pj ; 

problem  Min_Scenarios_N  :  z,  No_Scenarios  ,  Factor_Value  ; 
sol  ve  Min_Scenarios_N ; 

display  z ; 

display  {j  Fn  N}  sum{  i  hi  S[j]}  z[i]; 

# - # 

#  Integer  Linear  Programming  (outcome  j  occurs  is  set  to  Pj ) 

#  - - - - - -ft 

problem  Min_Scenarios_P  :  z,  No.Scenarios ,  Factor_P; 
solve  Min_Scenarios_P ; 
display  z ; 

display  { j  in  N}  sum! k  _m  S  [  j  ]  }  yave  [ k ]  *  z  [  k  ] ; 

#  -  - # 

#  Integer  Linear  Programming  (combining  both  —  not  illustrated) 

#  - - - - - - - - -  - # 

problem  Min_Scenarios_P_N :  z,  No_Scenarios ,  Factor_Value  ,  Factor_P ; 
solve  Min_Scenarios_P_N ; 

display  z ; 

_ / 
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